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Using the recently published data of twin kHz quasi-period oscillations (QPOs) in neutron star low-mass X-ray binaries 
(LMXBs), we study the different profiles between bright Z sources and less luminous Atoll sources in a statistical way. 
We find the quality factors of upper kHz QPOs show a narrow distribution both for Z sources and Atoll sources, which 
concentrate at 7.98 and 9.75 respectively, the quality factors of lower kHz QPOs show a narrow distribution for Z sources 
and a broader distribution for Atoll sources, which concentrate at 5.25 and 86.22 respectively. In order to investigate 
the relation between the quality factor and the peak frequency of kHz QPOs, we fit the data with power-law, linear and 
exponential functions, respectively. There is an obvious trend that the quality factors increase with the peak frequencies 
both for upper and lower QPOs. The implications of our results are discussed. 
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1 Introduction 

A number of neutron stars in low mass X-ray binaries show 
the kilohertz quasi-periodic oscillations in their X-ray spec- 
tra (van der Klis 2000, 2006). These frequencies, in the 
range of 200 ^ 1300 Hz, are as the same order as the dy- 
namical timescale of the innermost region of the accretion 
flow around the stellar mass compact objects (van der Klis 
2006, 2008). Owing to the expected links with the orbital 
motion, most work about kHz QPOs focus on the expla- 
nation for the nature of these signals (e.g. Miller, Lamb 
& Psaltis 1998; Stella & Vietri 1998, 1999; Kluzniak & 
Abramowicz 2001; Abramowicz et al. 2003; Zhang 2004). 
The kHz QPOs often arise as simultaneous twin peaks (up- 
per vi and lower v\ frequency) with frequencies chang- 
ing over time. These frequencies behave in a rather regular 
way and follow the tight correlations between their frequen- 
cies and other observed characteristic frequencies (see, e.g. 
Psaltis et al. 1998, 1999a; Psaltis, Belloni & van der Klis 
1999b; Stefla, Vietri & Morsink 1999; Befloni, Psaltis & 
van der Khs 2002; Titarchuk (& Wood 2002; Mendez (& van 
der Klis 1999, 2000; Mendez et al. 2001; Yu, van der Klis, 
Jonker 2001; Yu, van der Klis 2002). Moreover, the correla- 
tion between the upper frequency and the lower frequency 
across different sources can be roughly fitted by a power- 
law function (see e.g. Psaltis et al. 1998, 1999a; Zhang et 
al. 2006), and also by a linear model (see Belloni, Mendez 
& Homan 2005, 2007). 

The kHz QPOs in LMXBs are narrow features (peaks) 
in their power density spectra (PDS), whose profiles can 
be described by the Lorentzian function Py oc A{iwl\(y — 
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vq)^ + (u'/2)^] (j/Q is the peak frequency, w is the full width 
at half-maximum (FWHM), and Aq is the amplitude of the 
signal). The ratio of peak frequency to FWHM is the quality 
factor, 

Q^^. (1) 

w 

Therefore, the kHz QPOs in neutron star LMXBs can be 
characterized by three characteristic quantities, i.e. centroid 
frequency (i.e., peak frequency i/q), quality factor (Q = I'o / 
FWHM) and fractional root-mean-squared (rms). The qual- 
ity factor characterizes the coherence of the signal, while 
the rms represents a measure of the signal strength and it is 
proportional to the square root of the peak power contribu- 
tion to the PDS. Each kHz QPO corresponds to its quality 
factor (lower Q\ and upper quality factor) and centroid 
frequency (upper and lower v\ frequency). 

In the past several years, the large Rossi X-ray Tim- 
ing Explorer (RXTE) archive makes it possible to study 
this quantity systematically in several sources. Using the 
data from RXTE, Barret et al. (2005a) studied the source 
4U 1608-52 and revealed a positive correlation between the 
lower frequency QPOs and its quality factors, up to a max- 
imum of about Q 200. Motivated by this idea. Barret, 
Olive & Miller (2005b, 2006) studied the QPO properties 
and the dependency of its quality factor on the peak fre- 
quency in source 4U 1636-536. It shows that the quality 
factors for lower and upper kHz QPOs of 4U 1636-536 fol- 
low different tracks in a Q vs. v plot, i.e. the quality fac- 
tors of lower kHz increase with its peak frequency up to 
850 Hz (Q ^ 200) and drop precipitously to the highest 
detectable frequency ^ 920 Hz (Q ~ 50), while the qual- 
ity factors of upper kHz QPOs increase steadily all the way 
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to the highest detectable frequency. Moreover, the quality 
factors of lower QPOs is higher than those of upper ones 
(Barret, Ohve & Miller 2005b,c; 2006). The rough similar- 
ity was also extended to sources 4U 1735-44 and 4U 1728- 
34 (Barret, Olive & Miller 2006; Boutelier, Barret & Miller 
2009; Mendez 2006; Torok 2009). 

In this paper, we analyze the distributions of the qual- 
ity factors for upper and lower kHz QPOs in ten sources 
— five Atoll sources (namely 4U 1608-52, 4U 1636-53, 4U 
1728-34, 4U 1820-30 and 4U 1735-44) and five Z sources 
(namely Sco X-1, Cyg X-2, GX 17h-2, GX 5-1 and GX 
340H-0). In order to investigate the relation between the qual- 
ity factors and its centroid frequencies, we fit the data with 
the three functions , i.e. linear relation {Q = a + 6z^/1000), 
power-law relation (Q = a(z^/1000)'') and exponential re- 
lation (Q ^ a X exp{bh'/1000)), according to Q — i' tracks. 
We benefit from the existing studies and use the published 
data from the collection by Mendez (2006, and reference 
therein). In section 2, we analyze the data for these sources 
and execute fitting for the different source categories. Con- 
clusions and discussions are contained in section 3. 

2 Statistical Analysis for the Coherence of 
kHz QPOs 

In this part, we investigate the data statistically and choose 
linear (Q = a + 6z^/1000), power-law (Q = a(z^/1000)'') 
and exponential relation {Q = a x exp(bu/1000)) for fit- 
ting, where a and b are the undefined parameters. Firstly, 
we put the sources together and fit them with above three 
relations, then we divide them into Atoll and Z sources and 
execute the same fitting. 

2.1 Statistically Analysis for the Data 

Firstly, we study the distributions of upper and lower quality 
factors forZ and Atoll sources. Fig.[T]presents the distribu- 
tions of upper quality factors {Q2) for Z and Atoll sources. 
It is found that most Q2 of these two classes locate in the 
similar region (i.e. (52 = 2 - 18) as a whole. For the excep- 
tion, a few Q2 of Atoll sources present relatively larger val- 
ues (up to 40, see Fig.[T]for detail). Most Q2 for Z sources 
gather in a range from 2 to 10 which is larger than those 
of Atoll sources. For a quantitative knowledge, we calcu- 
late the mean value of Q2 for Z and Atoll sources, i.e. < 
Q2Atoii >= 9.75 and < Q2Z >= 7.98. The mean value of 
Q2 for Z sources is lower than that for Atoll sources. 

In Fig. |2] we plot the distributions of lower quality fac- 
tors (Qi) for Z and Atoll sources, where we notice very big 
difference in the range of Qi between these two sources. 
The Qi for Z sources are small {Qiz < 14), and those 
for Atoll sources are very large (up to Qi = 200). The Qi 
ranges from 2 to 14 for Z sources and from 2 to 200 for Atoll 
sources (see Fig.|2]for detail). Most data of Qi for Z sources 
are distributed from 6 to 10, and that of Atoll sources gather 
in the region of Qi ^ 60 — 100. Some Qi for Atoll sources 
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Fig. 1 The distributions of Q2 for Z and Atoll sources. 
The red shadow represents the Q2 distribution for Z sources, 
while the black shadow show that for Atoll sources. The 
mean value of Q2 is 7.98 for Atoll sources and 9.75 for Z 
sources. For the plot, we take the width of bin as 2. 

are as high as 200. We also calculate the mean values of 
Qi both forZ and Atoll sources, i.e. < QiAtoii >~ 86.22, 
< Qiz >= 5.25. 

On the whole, the quality factors for Z sources are smaller 
than those for Atoll sources. In addition, Q2 for Z sources 
cover a wider range than Qi, while Q2 cover the narrower 
range than Qi for Atoll sources. The mean values of Qi and 
(32 are < Qi >= 49.20 and < Q2 >= 8.70 , respectively. 

2.2 Fitting for the Q2 1^2 Relations 

In order to study the Q2 — V2 relation , we use linear (Q = 
a + 61//IOOO), power-law {Q = a{v/lQQQf) and exponen- 
tial relations (Q = a x exp{bv / IQQQ)) to fit Q2 — V2 tra- 
jectories for ten sources. It is noticed that the Q2 values are 
similar for Z and Atoll sources, which increase with V2, so 
we fit the three relations to total ten sources, Z sources and 
Atoll sources, respectively (see Fig. [3]). The values of fitting 
parameters (a and b) and the fitting results are listed in table 
ffl 

From Fig. [3] it is seen that the data points are dispersive 
in the left and right panels, while the result in the middle 
panel is better, especially for power-law and exponential re- 
lations. From table[T] we notice that the fitting results are not 
good for these large /d.o.f.. If we just forecast a roughly 
trend based on these results, it can be seen that the power- 
law and exponential relations fit mildly better than that of 
linear relation as a whole. Exponential relation fits better 
than power-law and linear relations both for Z sources and 
for Atoll sources. The fittings for Q2 — 1^2 of Z sources are 
better than that for Atoll sources, with relatively small re- 
duced ix^ /d-O-f- = 4.73, 5.01 and 5.76 for exponential, 
power-law and linear relations, respectively, see table [T] for 
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Fig. 2 The distributions of Qi for Z and Atoll sources. The left panel is for Atoll sources, and the right one is for Z 
sources. The mean value of Qi is 86.22 for Atoll sources and 5.25 for Z sources. For the plot, the width of bins are set to 
20 and 2 for Atoll and Z sources, respectively. 



detail). But for five Atoll sources, the reduced are rel- 
atively large (x^/d.o.f. ^ 24.63, 26.10 and 29.33 for ex- 
ponential, power-law and linear relation, respectively), and 
the exponential fitting is also better than the other two re- 
lations. If we put all the sources together and fit three re- 
lations to them, we find that the power-law relation with a 
small x^/d-o.f. (1 1.75) is better than the other two. 

2.3 Fitting for the Qi - Relations 

Following the same techniques, we investigate the Qi — i^i 
relations for Z and Atoll sources (see Fig.lU. 

The left panel of Fig. |4] show the results for Z sources, 
in which it can be seen that almost all points form the rising 
trend, and we fit power-law, linear and exponential relations 
to all the data points, respectively. The parameters and fit- 
ting results are listed in table |2] The reduced ^re 3.33, 
3.30 and 3.62 for power- law, linear and exponential rela- 
tion, respectively. It seems that the linear relation fits better 
than the other two. 

The right panel of Fig.|4]show the results for Atoll sources, 
in which it can be seen that the Qi — vi distribution of 
five Atoll sources display two rising parts and two drop- 
ping parts (see the right panel of Fig. |4]for detail). We take 
the sources 4U 1608-52 and 4U 1820-30, which have large 
Qi values, as a group. While the sources 4U 1636-53, 4U 
1728-34 and 4U 1735-44, which have small Qi values, as 
the other group. Then we fit the data of the two groups re- 
spectively, the fitting results are in table |2] Compared with 
the fitting result of Q2 — 1^2, the reduced is very large 
(see table |2] for detail), but also with large R^. It is found 
that the linear relation fits better than the other two for the 
curve with relatively low Qi values (formed by 4U 1636- 
53, 4U 1728-34 and 4U 1735-44), and the power-law rela- 



tion fits better to the curve with high Qi values (formed by 
4U 1608-52 and 4U 1820-30). 

So far, there is no appropriate model which can explain 
the mechanism of lower kHz QPOs, as well as that of the 
lower quality factors. However, it is claimed that the lower 
kHz QPOs arise in the innermost region of accretion disk 
(van der Klis 2006). The quality factors relate to the fre- 
quency drift in the inner disk (Barret, GUve & Miller 2006; 
Wang et al. 201 1). So we claim that the diverge of the data 
points between different sources may result from different 
physical environment of inner disk in different sources. In 
addition, the occurrence of drops in Qi — i^i plots for Atoll 
sources may arise from the existence of the inner bound- 
ary of accretion disk (Barret, Olive & Miller 2006). Due to 
the different mechanism of the drops of Qi — 1^1 tracks for 
Atoll sources, we just fit the rising branches of the Qi — ly'i 
track and do not consider the dropping parts with the same 
techniques. 

3 Discussions and Summary 

In this work, we study the distribution of the quality fac- 
tors and the relation between the quality factor and peak 
frequency of kHz QPOs . We notice that there are differ- 
ences in the distribution of the quality factors between lower 
and upper kHz QPOs, as well as Z and Atoll sources. By fit- 
ting the power-law, linear and exponential relation to the 
tracks, we investigate the relations for upper and lower 
kHz QPOs. 

(1). The Q2 values are low {Q2 < 18 in general, with a 
few exceptions) for both Z and Atoll sources. The Qi values 
for five Z sources are very low, Qi^ = 2 — 14. However, 
the Qi values for five Atoll sources are very high, and the 
maximum QiAtoii is up to 200 (4U 1608-52). 
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Fig. 3 The fitting plots of Q vs. v for upper kHz QPOs. The left panel is for all ten sources. The middle (right) panel 
is for five Z (Atoll) sources. The three parts from top to bottom of each panel present the power, linear and exponential 
fitting, respectively. The point with different shapes and colors in the plots represent the different sources as shown in the 
legends. The fitting functions and results are shown in table[T] 

Table 1 The fitting results for quality factors of upper kHz QPOs correspond to Fig. 3 . 





Function 


a 


b 








Total fitting 












= a + bu2/W00 


-9 ±2 


21 ± 2 


14.00 


0.47 


Q2 


= a(i/2/1000)' 


10.9 ±0.3 


2.2 ±0.2 


11.75 


0.49 


Q2 = 


a X exp{bu2/1000) 


0.9 ±0.2 


2.5 ±0.2 


11.80 


0.48 




Z source 










Q2 


= a + bu2/W00 


-10 ± 1 


21 ± 2 


5.76 


0.64 


Q2 


= a(i^2/1000)* 


12 ±0.4 


2.8 ±0.3 


5.01 


0.69 


Q2 = 


a X exp{bu2/1000) 


0.4 ±0.1 


3.4 ±0.3 


4.73 


0.71 




Atoll source 










Q2 


= a + 6!/2/1000 


-13 ±5 


25 ±5 


29.33 


0.36 


Q2 


= a(i^2/1000)' 


10.9 ±0.9 


4 ±0.8 


26.10 


0.43 


Q2 = 


a X exp{bu2/W00) 


0.1 ±0.1 


4.3 ±0.8 


24.63 


0.46 



(2) . On the whole, the ranges of Q for five Atoll sources 
are wider than that for five Z sources both for upper and 
lower Q. According to the idea that the twin kHz QPOs 
come from the inner region of accretion disk, the instabil- 
ities which can lead to the radial drift and frequency drift 
in this region for Atoll sources are stronger than that for Z 
sources (Wang et al. 2011). 

(3) . From Fig. |2] it seems that the range of Qi for Atoll 
sources is about 10 times wider than that for Z sources. In 
addition, the high luminosity Z sources exhibit luminosity 



close to critical Eddington luminosity L Edd, and the range 
of their luminosity is 0.5 — IL Edd- But Atoll sources present 
low luminosity 0.001 — Q.2LEdd, 100 times lower than that 
for Z sources (Hasinger and van der Klis 1989; Hasingger 
1990, see van der Klis 2006 for a review). Therefore, we 
expect a relation of Qi ~ L"^/^ between the luminosity 
and quality factors for lower kHz QPOs to be investigated 
and proven in future. 

(4). The fitting results of Q2 vs.V2 can not confirm the 
actual relation between Q2 and V2 because of the large /d.o.f.. 
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Fig. 4 The similar pattern as Fig. |3] but for lower kHz QPOs. The left panel is for five Z sources. The right panel is for 
five Atoll sources, in which we take the sources 4U 1608-52 and 4U 1820-30 (higher quality factors) as a group and the 
sources 4U 1636-53, 4U1728-34 and 4U 1735-44 (lower quality factors) as the other group, then we fit the data of two 
groups, respectively. The fitting functions and results are shown in Table. |2l 

Table 2 The fitting results for quality factors of lower kHz QPOs correspond to Fig.4. 



Function 


a 


b 


x'/DoF 


R' 


Z source 










Qi = a + bui/WOO 


-1.3 ±0.6 


12 ± 1 


3.30 


0.59 


Qi = a{ui/WOO)'' 


11.3 ±0.9 


1.3 ±0.2 


3.33 


0.59 


Qi=ax exp{bui/1000) 


1.4 ±0.3 


2.4 ±0.3 


3.62 


0.56 


Atoll source with lower Qi 










Qi = a + bui/lOOO 


-300 ± 30 


500 ± 30 


241.70 


0.87 


Qi = a{ui/1000)'' 


260 ± 30 


4.4 ±0.4 


343.75 


0.82 


Qi = ax exp{bvi/lOQO) 


1.1 ±0.5 


5.6 ±0.6 


380.99 


0.80 


Atoll source with higher Q\ 










Qi = a + bui/lOOO 


-200 ± 30 


480 ± 40 


270.52 


0.85 


Qi = 0(1/1/1000)* 


330 ± 30 


2.8 ±0.2 


267.59 


0.85 


Qi = a X exp{biyi/1000) 


7±2 


4 ±0.3 


279.81 


0.84 
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We consider that it may result from the small sample with 
the relative big errors. There is an obvious trend that Q2 in- 
creases with V2 for both Z and Atoll sources, and the power- 
law and exponential relations fit better than the linear rela- 
tion. 

(5) . The fitting results of Qi vs.vi also cannot be con- 
sidered as perfect. Then, an obvious trend is that the Qi 
increase with vi in the case of five Z sources, and the linear 
relation fits better than the other two. 

(6) . In the Qi — vi diagrams for Atoll sources, all five 
sources present drops at a maximum Qi. It is considered as 
the results of innermost boundary (ISCO) of accretion disk, 
firstly proposed by Barret and coauthors (see. Barret et al. 
2005b, 2005c, 2006, 2007) . The special mechanism of the 
drops are expected to be investigated further. Here, we just 
fit the three relations to the rising parts.) 

(7) . As for the relation between Qi and vi in Atoll sources, 
it can be seen in the right panel of Fig.|4]that, if only the ris- 
ing part of Qi — vi plot is considered, the five Atoll sources 
form two curve-like tracks, but it is not enough obvious to 
indicate a clear correlation. We claim that this phenomenon 
arises from the special physics in the inner disk region for 
each source, where stellar mass, magnetic field and accre- 
tion rate may be the roles of producing them. For the fit- 
ting, we fit the three relations to each rising curve. This 
tentative fitting is not good enough to show a clear corre- 
lation although the /d.o.f. of linear fitting is better than 
the other two for the curve with relatively low Qi values 
(formed by 4U 1636-53, 4U 1728-34 and 4U 1735-44), and 
the /d.o.f. of the power-law fitting is better to the curve 
with high Qi values (formed by 4U 1608-52 and 4U 1820- 
30). 
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